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SUW4ARY 

Weighti;ig is a well known technique for shaping the compresst.i pulse wavefonr 
it. t-adar processing. Usually, weighting is applied to the transfer I'unction 
and has the effect of sacrificing the time nainlohe width (resolution) in 
exchange f 'r decreasing the height of the neighboring siielobes. This paper 
reports on simulations of weighting in the time dom*'lii, as used to shape the 
lime-compressed pvJ.se waveform, nie digital input radar data is 32 bit 
I,Q, and simulates data from a point target as imaged by a Seasat-A type 
system. Weighting functions tested include stepped-amplitude distributions, 
(with 1 through 5 steps), anc the cosine-squared plus pedestal distribution. 
Effects treated include mainlobe broadening, peak energy reduction, the 
intigrated sidelobe rp.tio, signal to noise ratio, and nearest sidelobe 
suppression. 

1.0 INTRODUCTION 

Pulse compression in radar is commonly accomplished through matched filtering. 

A radar point target return s(t) is correlated against a T-shifted version 
of its conjugate s*(t) over a time interval T, to yield an output at tiune t: 

rT/2 

s(x) =1 s(t)s*(t - t) dt. (l) 

J_T/2 


• This paper presents the results of one phase of research carried out at the 
Jet Propulsion Laboratory, California Institute of TechnoD.ogy, under Contract 
No. NAS7-100, sponsored by the National Aeronautics and Space Administration. 
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The estimate s(t) yields a prediction of the return from a target at the 
location specified by t. The function 1s(t)| usually has multiple peaks, 
and thus creates detection ambiguities when there is la7'ge signal dynamic 
range. For example, if s(t) is a linear P^^ signal, s(v) approximates e sine 
function with high 1-13.2 db) sidelobes. 

The relation (l) may also be expressed in the frequency domain. Set 

s (t) = s*(-t). 
r 

Then (l) becomes 


rT/2 

§(t) = I s(t) s^(t - t) dt (2) 

•'-T/2 

= s0s^(t). 

If capital letters denote fourier transforms, then (2) Implies 

S(w' = S(w) • Sp(v). (3) 

The inverse transform of F(w) is of course s(t). In (3) we can multiply 

S (w) by a weighting function W(w) to obtain a new output S_ : 
r 1 

S. (w) = S(w) • S (w) •*W(w). (1*) 

i r 

The inverse transform, s^(r)» is now only an approximation to s(t). The goal cf 
frequency domain weighting is to select a function W in (l») such that s^(t) is 
a good local approximation to §(t), but has lower sidelobes. 

It is known that weighting of the time domain signal can produce similar 
effects, particularly if the signal iu linear FM. In time weighting the 
weighting function is multiplied with the return signal 
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s^(t) = s(t) wft) 

or with the reference function s"*(t) 

3*(t) = s*(t) w(t), 
prior to filtering. 

■fliis paper reports on simulations time-weighting, performed on digital 
SAR data of the type expected from satellites such as SEASAT-A. 

It is showii that use of a cosir.e-squared-plus-pedestal time-domain weighting 
scheme can result in sideiobe suppression down to -36 db, with mainlobe 
broadening of less than 30 % for the simulated data. 


Section 2 gives a brief description of the frequency domain weighting approach, 
which is th«;n used to motivate the ideas on time weighting presented in 
section 3. Finally, section 3 contp.ins the numerical results extrapolated 
from the simulations. 

2.0 FREQUENCY DOMAirj \7EIGHTING 

Let the signal given by eq. (5) be received at the radar receiver at time t. 


s(t) = 



The matched filter response at time t is 


( 5 ), 



-j/|(t+T )^+f^(t+X )j 


dt . 


( 6 ) 
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By elementary calculation, we have 


s(t) = 2Te 


-^v 


sin (bTT - 


bTi 


(7) 


Usually, the range of interest is such that t is small compared to T, so that 
the approximation 


s(t) 


2Te 


-JV 


sin (bTi ) 
bTt 


( 8 ) 


= 2T sine (bTx) x (phase factor) 

is valid. Ignoring the phase factor, we see that s(x) peaks at x = 0. There- 
after, sidelobes 13.2 db below the peak intensity occur near x = ± 3^/2bT. 


These sidelobes can mask weaker peaks corresponding to neighboring targets, 
as shown in Figure 1. 

By slightly mismatching the filter response, (and hence sacrificing signal to 
noise ratio and nainlobe width), the sidelobe levels of the filter output can 
be decreased. The process is most easily itnderstood when mismatching is 


-«»-T 



Fig. 1, Sidelobes Masking Weak Peaks 
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. \ 

vHei'e B is proportional to T, ny virtue of (8)). When the three sine func cions 

are nlded, a new "quasi-sinc" function results with lower sidelobes. (tJee fig-- ’ 

ure c). 

Thus weighting by W(-rf) = [1 + cos 2nw/B] can lower the time sidelobes. In 
the radar case, we would actually be confronted with a superposition cf sine 
functions at the vuiweighted filter output. Each sine function would bo centered 
on its own resolution element. However, by the linearity of the fourier trans- 
form operation, one weighting function applied to the superposition s)>ectrum 
serves to suppress the sidelobes of all sine functions u!-, once. 








3 . 0 TIME DOMAIN WEIGHTING FOR TIME SIDE:j03E SUPPRESSION 

Time domain weighting is accomplished by multiplying the return signal or time 
reference function by a propei’ly designed -weighting function. Since the weights 
can easily be incorporated as part of the time domain reference function, the 
technique is a promising one for application to real-time spacecraft on board 
processors, 'the effect of time weighting on time output can often be easily 
determined [l]. 

For example, let the radar signal s(t) be the linear FM signal given in (5). 

Then, in a rough manner of speaking, there is a one-to-one correspondence between 
time and frequency. That is, each time interval At corresponds to a frequency 
interval Aw in the spectrum. Hence a time weighting function W(t) would have 
about the same effect as the frequency weighting function W(t = w). Thus a first 
order estimate of the effects of time weighting by W(t) may be obtained through 











frequency veifihting with W(w). More precisely, it is known that for the 
linear FM signal (5)» the filter time output g(t) resulting from a tine 
weighting w(t) , satisfies 


ls(t)l 



w(t) exp {JbtT} dr. 


For further details, see [l] p. 18Y. 


(13) 


In SAR processing, the azimuth aperture is generated by a sensor moving linearly 
with constant velocity over the upertui'e length. Thus, to u reasonable 
ajjproximation, the aziiuuth matched filter function Isa linear KM chirp, and the 
above time weighting approximations (e.g. (13)) aoply. However, tlie azimuth chirp 
is a sampled (discrete) function, and tlie results ai e somewhat different. 
Therefore, to discover the effects of weighting with discrete time functions, 
actual simulations were perfomed. The v.'eigiit ings wore applied along tiio azinutli 
direction of simulated f'RASAT-A SAK data. Kach data inpr.b point wns quantized 
to bits, us was each reference and weighting coefficient. Tiic gro\ind spacing 
between adjacent data input points was around meters. A total of -’iCOC input 
points was used to process a b-look output data po'nt. l.hus each look required 
filtering of 1024 points, or 1024 separate weighting coefficients. The 
spacing between output data points was around l6 meters. The simi.ia'..ed input 
data was designed to represent the retui-n from a point target so that the effects 
of weigliting could be easily observed. The output image was ]6 (range) by 512 
(azimuth) points. 


The weighting functions selected for simulation vfore i.he stepped :implitude 
distributions and the co.sine - squared plus pedestal function [3], [4], [5]. 

The stepped amplitude distributions are given in table 1 for 1 th:'ough 5 steps. 
(Gee also figure 3). These distributions were chosen for their optimality in 
antenna pattern adjustment, and their ease of implementation (3]. sote that 
T is proportional to N = 1024, the number of weight.^-. 
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TABLE 1 


STEPPED Al-IPLITUDE TIME WEIGHTING FUNCTIONS 



^2 



a b 

5 1 

^2 

■''3 



1 

0.5 

0.5 

% • * 

!!! 

1 

1 

0.55 



• • • 

0.35 

0.35 

0.30 

. . . 

1 

O.b'25 

0.350 

. » , 

. • • 

0.25 

0.25 

C.25 

0.25 

1 

0.78 

0.56 

0.3»i 

• 

0.300 

0.225 

0.235 

0.170 

0.070 1 

0.72 

0.5^< 

0. S6 

0.18 


(— 



Fig. 3- Stepped Amplitude Time Weighting for Time Waveform Shaping 

The cosine squared-plus-pedestal distribution [5] ic, through a tri>. ..ometric 
identity, equivalent to the weighting given by the right-hand factor in eq. (lO). 
It is defined by 


W(t) = H + {1 - H) cos^ (irt/T) 


(14) 
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where 


H = pedestal height (O ^ H < l) 

T = pulse timewidth. 

For the simulation, ve select H = O.OS to approximate Hamming weighting, which 
produces the lowest sidelobes attainable ./ith this type weighting in the 
frequency domain [5]. The resulting weighting is similar to the Taylor 
approximation to the physically unrealizable optimum Dolph - Chebyshev weighting 
[2], [6], and represents a practical approach for digital processing. 

Table L gives the total energy in decibels in the 3 center lines about the peak. 
(By "line" is meant a range line of l6 data points centered around the azimuth 
peak. See figure 


TABLE 2 

PEAK ENERGY AND INTEGRATED SIDELOBE FJ^TIO 


Weighting Technique 

Total Energy in 3 Center Lines 

Peak Energy 

Rectangular 

103.109 

100.067 

2 step amplitude 

103.»*58 

99.66? 

3 step amplitude 

103.51^ 

99.1*15 

I4 step ajnplitude 

103. 5 Vc 

99.359 

5 step amplitude 

103.552 

99.362 

003^*.+ pedestal 

103. >*78 

98.71*8 


Note: Total energy in 512 lines is 103. 671 db for each type of weighting. 
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O O O O O 

PEAK 

O O O ® O O 

O O O O O 

RANGE »- 

Fis. 

Since the total energy in each l6 x 512 image was normalized to 103.671 db, 
the energy in the center lines reveals in a rough way how much energy is left 
over for sidelobe generation. This data can of course be used to compute the 
"integrated sidelobe ratio" in decibels, defined by 

total enei'gy outside mainlobe 
total energy in main lobe 
(3 center lines) 



Integrated Sidelobe Ratio = 10 log ^ 

J t) 


® . AZIMUTH 

O 

O 


The results 
give a fair 


Rectangular 

(no weighting) 

2 - step amplitude 

3 - step amplitude 
h - step amplitude 
5 ~ step amplitude 
cos^ + pedestal 


for both U and 32 bit quantization, and 


DB FOR V.^RIOUS FIGHTINGS 


Integrated 

Sidelobe Ratio 

b bits 

32 bits 

- 12.1 db 

- 12.8 db 

- 13.3 db 

- IU.3 db 

- ll*.5 db 

- 16.1 db 

- 15.6 db 

- 17.8 db 

- 15.7 db 

- 18.0 db 

- 17.3 db 

- 21.7 db 


are summarized in table 3, 


estimate oi’ the 2-dimensional integrated sidelobe ratio. 


TABLE 3 

INTEGRATED SIDELOBE RATIO 


Weighting 
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l':ib.lc C n.l;;o an itKiioation o’’ loss in siKiial to noise »*atio in icras of 

peak encr,;i'. iUnce eai’h weighted output has the same total onorj'y, the peak 
eiuM>',y in the second Cv'lunm indicates the loss in PNK rather aoeuratoly. The 
^'.t•eatest loss (l.3l9 db) reiiitive to the unweij^hted case occurs vitli cos' + 
pedoi tal veljr.htiiu:. 

Table u is used to derive the inainl.v'>be breadeninp resulting, fro:;i vei|r.hti;ir.. As 

% 

o\pectod, the v.orst broadeniuf; is seen to occur with cos ' + ped vei»:htiii(',» w'l;l’ 
only a 1 dh suppiams.ion l6 meters from the peak. Tii this case the i .ib ro olutic!i 
is .lc;:raded to oO meters. 

Tb.e values in Table !i ean be \ised to derive faiviy preoiso values for the 
ri.'solutiou achieved witli each wv iplitinp scheme, in the folloviiift manner. The 
.j'.it.put. filler fi.hctiv'ii is tiiovie’JeJ as a sine Ameti(.’n wliich is clcsely ti’ack'd 
in i.ho ratjpe of interest by a quadratic. Usin>: Laf^ran.ce’ s interpolation pi iy- 
nom'al, with the peak and .1 nearest nei*thbors as input- points, the distance to 
t.h.e 0 ub mainlv-d'o threshold is ealeuiated. The I’csxilLs are shown in table ?. 

’FAhliK 

at fjiiVKKfT kiadinvH, pi? 


'.\'c if.l’.t i ap. T’eciini quo 


ttuppressiv-'u , db 


Hect.'U'.ftular 

step a.mpli tilde 
i step -uiipli tilde 
•’i step loplitndo 
'} step ainplitude 
Cvis” pedestal. 
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TABLE 5 


RESOLUTION ACHIEVED WITH WEIGHTING 


Weighting 

Resolution 

Rectangular 


(No weighting) 

?.l m 

2 - step :unplitude 

27.2 m 

3 - step amplitude 

28.8 m 

•» - step.:implitude 

28. 8 m 

5 - step a’lplituJe 

28.8 ra 

Cos‘" + pedestal 

30. L m 


Kinally, table (i gives the siclelobe suppi'ossion actually achieved with each 

weighting, scheme. Each technique generated 2 side.lob?s of generally different 

heights; hence the min imd max. suppressions are both tabulated. The best 

o 

sr.ppress ion is ne.liieved by cos" + pedestal weigliting, with a -3o db measure- 
ment at the highest sidelobe. (The sidelobe asjsnmetry results from the 
true peak v>ceurring berween two data output poii\t.s). 


TABLE 6 

NE.AKKST SIDELOBE SUPPRESSION, DB 


Weighting Toch:iique 

Sidelodc Suppreti 

ision, ' 


min 

Max 

Rectangular 

20 

-2h 

2 step amplitude 

22 

2S 

3 step ampl itude 

26 

29 

L step amplitude 

30 

33 

5 stop amplitude 

29 


Cos*' + pedestal 

3b 

!*0 
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5i,o cor;cLi]:Mo:i 

Ve tine k;.\T.;iin I’or tine si Jejooe suppression pivsent.i' ;i:i attraj- 
t ivo alternative to Trequenoy aonuin weiyju iiit';, espeeial ly if tijne-.iona in 
auitnali; processing :s uo be ucconplishoii . 'i'tie simuiationr. show tiint cos’ + 
iva-o.jtui weighting aoi.ieves 3o db sLJelobe suppression with a GNU joss of less, 
t.'.'.au i.5.vib, and resolution loss of around ?5o. 


Tae nnthor would like to th:mk Dr. f. Ku for valuable software a:ui aisoussions , 
and also H. i’iereson, W. Arens, and V. Tyi'ec for helpful suggest! tms . Irv.nlu.-iL 
progra'.nr.ing aid was furnisheri by B. Barkan. 
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